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coefficient in the range of energies where transmission resonances occur. We calculate 
the Wigner time-delay and show how their maxima depend on the position of the 
transmission resonance. 
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1. Introduction 

Barrier penetration of relativistic electrons is a very important problem in scattering 
theory and provides a theoretical framework for different physical phenomena that 
are not present in the non relativistic regime such as the Klein-paradox and 
super criticalityjU EJ [3] 

The discussion of tunneling of relativistic particles by one- dimensional potentials 
has been restricted to some simple configurations such as delta potentials and square 
barriers, mainly in the study of the possible relativistic corrections to mesoscopic 
conduction jl] and the analysis of resonant tunneling through multi-barrier systems [5]. 
Recently[6], electron transport through electrostatic barriers in single and bi-layer 
graphene has been studied using the Dirac equation and barrier penetration effects 
analogous to the Klein paradox. 

The study of transmission resonances in relativistic wave equations in external 
potentials has been extensively discussed in the literaturej2j [7J El [9]. In this case we 
have that, for given values of the energy and of the shape of the barrier, the probability 
of transmission reaches unity even if the potential strength is larger than the energy 
of the particle, a phenomenon which is not present in the non relativistic case. The 
relation between low momentum resonances and supercriticality has been established 
by Dombey et al^\10\. Recently [9| fTU| fTTj . some results on scattering of Dirac particles 
by one a dimensional potential exhibiting resonant behavior have been reported. 

The study of the tunneling effects of Dirac particles by potential barriers has been 
almost restricted to those cases where the wave equations are solvable in terms of special 
functions and a straightforward identification of the asymptotic states is possible. The 
composition of barrier potentials [12] using the scattering matrix method permits us to 
consider physical configurations which cannot be solved in closed form. We discuss 
resonant tunneling of Dirac particles by a double square barrier and double cusp 
potential when transmission resonances are also present. We show that the presence 
of transmission resonances close to the position of the poles of the scattering matrix 
turns out in the appearance of additional peaks in the transmission coefficient and in a 
modification of the Breit-Wigner Lorentzian profile. We calculate the Wigner time-delay 
and show that some of the peaks cannot be associated with energy resonances. 

The article is structured as follows. In Sec. 2, we calculate the transmission of a 
Dirac particle by a double barrier of equal strength in a range of energies where the 
barriers support transmission resonances. We calculate the Wigner time-delay of the 
resonances and show that, the transmission coefficient as a function of energy exhibits 
peaks that cannot be identified as energy resonances although they correspond to peaks 
of the Wigner time-delay. In section 3, applying the composition of scattering matrices, 
we discuss the transmission of a Dirac particle by a double cusp potential. We show 
that when the system exhibits energy resonances in the range of energy where the cusps 
support transmission resonances the transmission coefficient and the Wigner-time delay 
also present maxima that cannot be identified as Breit-Wigner resonances. In Sec. 4 
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we briefly summarize our results. 



2. Double square barrier 

Resonant scattering of a Dirac particle by a square potential barrier has been discussed 
by different authors in the literature [3]. Here we are interested in studying resonant 
transmission of relativistic particles by two square barriers separated a distance d when 
barriers support transmission resonances [lj. 

Since we are working in (1 + 1) dimensions, we choose the following representation 
for the 7 matrices: 

where o~\ and 03 are the Pauli matrices. Using the representation ([T]), the Dirac equation 
in the presence of a potential V(x) takes the form[TT] 

[a 2 (E - eV(x)) + a l d x + m]tf = (2) 

where we have adopted the natural units h = 1 , c = 1 . 

The solution to the Dirac equation in the presence of a potential barrier of height 
V can be obtained after decomposing the spinor in three regions. For x < (Region I), 
the solution to the Dirac equation ([2]) has the form 

*z = A 1 ( l_ ) e fe + B, ( X_ ) (3) 

\ E-m ) \ E-m J 

The spinor solution in the sector II (0 < x < a) is 

*!! = <*[ I ) e ipx + P [ X \e*> (4) 

E-m-eV I \ E-m-eV 



where p = ((E — eV) 2 — m 2 ) 1 / 2 . In the region III x > a we have 

*// = Ci ( _i_ ) e ik » + D X ( _X ) e~ lkx (5) 

\ E-m J \ E-m J 

Imposing the continuity of the spinor solution at x = and x = a, we obtain that the 
relation between the amplitudes of the incoming and outgoing waves is 

B = -2i(l - 7 2 )sinpq ^ ^e~ lka ^ 

1 (1 - 7) 2 e^ a - (1 + 7) 2 e - i P a 1 (1 - 7) 2 )e ipa - (1 + 7) 2 e" i P a 1 

4je~ ika ^ -2i(l - 7 2 ) sinpae~ 2ika ^ 

1 ~ (1 - 7) 2 )e^ a - (1 + 7) 2 e"^ a 1 (1 - ^) 2 e^ a - (1 + 7) 2 e"^ a 1 

Using Eq. (jBJ and Eq. (CO) we readily obtain that the scattering matrix S\ associated 
with a square barrier of strength V in [0, a] is 

" 2 ) sin pa _4^ e - lfea 

-ika OoYl ^/2"\ n^nr,— 2ika 



11 — 2i(l — 7 2 ) sin pa — 47c" 

bi - 



(1 - 7) 2 e'P a - (1 + 7)26-*° \ -4 7 e" iA:a -2i(l - 7 2 ) sinpae" 



(8) 
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where 7 = ^-pp^^ • Applying the property of transformation of the scattering matrix 



under a translation [12J, we obtain that the matrix S 2 associated with a potential barrier 
of strength V\ in [d, d + b] is 

' -2ie 2ikd (l - 5 2 ) sin qb -45e~ ikb 



S 



§y e iqb — (1 + 5) 2 e~ iqb \ -A5e~ ikb -2ie~ 2ikd (1 - 5 2 ) sin qbe~ 2ikb 



(9) 



where 5 = ^^E-m)^ ■ The composition of two scattering matrices Si and S2 can be 
written as|13] 

s = < = «r ,ri /+s& (10) 

where the indices 1 and 2 correspond to the potentials V and V\ respectively. 

Composing the matrices S\ and S 2 we obtain that the relativistic double barrier 
exhibits energy resonances for the values of E satisfying the equation 

((1 - 5) 2 e iqb - (1 + 5) 2 e~ iqb )({l - 7 )V pa - 

-(1 + 7 )V ipa ) + 4(1 - 7 2 ) sinpa(l - 5 2 ) smqbe 2tk(d - a) =0 (11) 

Transmission resonances for the double barrier occur when the transmission 
coefficient T of the whole system is equal to unity, this condition takes place, for a = b, 
when p = q and sin(pa) = 0, i.e. when the two barriers have equal strength V — V\ and 
the energy E satisfies the relation: 

E — V — yJn 2 -K 2 /a 2 + m 2 (12) 

Transmission resonances of a single barrier are maxima of the transmission coefficient, 
they are are not poles of the scattering matrix and therefore cannot be associated with 
quasibound states of the system, nonetheless their presence modifies the profile and 
peak distribution of the transmission amplitude against the energy. 

The relativistic double barrier exhibits a discrete number of resonances, whose 
position and shape reduce to those obtained in the double delta configuration as a — > 
and 6—7-0 with aV — > fi and bV\ — > 5 where /1 and 5 are the strength of the delta 
barriers. Since the transmission resonances for a delta barrier do not depend on the 
energy, all maxima of the transmission coefficient for the double delta are associated 
with energy resonances [T4"l [T5]. 

Fig. [2] shows the behavior of the transmission coefficient T versus the energy for 
two square barriers with equal width a = b = 3 and separation d = 5. The strength 
of the barriers has been chosen to allow the presence of transmission resonances across 
the barriers according to Eq. f fl2|) and the mass m has been equated to unity. The 
solid line in Fig. [2] depicts the case with V = V\ = 5 and the dashed line corresponds 
to V = 5 and V\ = 4. The poles of the scattering matrix S for V = V\ = 5 for 
1 < U(E) < 3 are E 1 = 1.5058 - 0.0688/, E 2 = 2.0414 - 0.0810/, E 3 = 2.7595 - 0.1141/ 
and the transmission resonances satisfying Eq. f fT2|) are located at E a = 1.7030 and 
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E b = 2.6791. Fig. [T] depicts the double barrier configuration with V = V\ = 5, d = 5 
and a = b = 3. It is also shown the position of the resonances for $tE < 3. From Fig. 
[2] it can be observed that there is a peak at the transmission resonance value E a . The 
transmission resonance at Eb is located close to E 3 and the peak does not exhibit a 
Lorentzian Breit-Wigner shape (see Fig. EJ). The dashed line corresponds to the case 
V — 5 and V\ = 4. It is worth mentioning the difference we can observe between the 
two plots around E b . 
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Figure 1. Double barrier configuration with V = V\ =5, width a = b = 3 and 
separated a distance d — 5. a) The dot line corresponds to dlEi — 1.505 b) The 
dashed line corresponds to the transmission resonance E a = 1.7033 c) The dash dot line 
corresponds to Sfti^ = 2.0414 d) The long dashed line corresponds to the transmission 
resonance Eb — 2.6791 e) The space dot line corresponds to 5Ri?3 = 2.7595. Notice 
that Eb is close to 5fti?3 



The phase shift (p of the transmitted amplitude t = \t\exp(iip) changes rapidly near 
the energy resonances. The Wigner time delay [16| [17] 

dip 
dE 



T 



(13) 
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Figure 2. Coefficient of transmission versus the energy for the double barrier 
relativistic system, with a = b = 3 and d = 5. The solid line corresponds to V = 5 
and Vi = 5. The dashed line corresponds to V = 5 and Vi = 4 

has maxima at positions very close to the resonance energies. The particle is trapped 
for a long time in region between the barriers before being transmitted. 

Figure 3 shows the Wigner time-delay corresponding to the resonances depicted in 
Fig. [2] The maxima of the Wigner time-delay are very close to the values of Ex, E 2 and 
E 3 There are also two additional peaks at E a and Eb that cannot be associated with 
energy quasibound states but exhibiting a behavior analogous to that one observed for 
energy resonances, resonances. Fig. [2] shows that the height the peaks at E a and E^ is 
smaller than those of the energy resonances between 1 < Re(E) < 3. 



Tunneling and transmission resonances of a Dirac particle by a double barrier 



7 




Figure 3. Time delay for the relativistic double barrier with a = b = 3 and d = 5. 
The solid line corresponds to V = 5 and V\ = 5. The dashed line corresponds to V — 5 
and Vi = 4. 



3. Scattering by a double cusp barrier 

In this section we proceed to study the resonant tunneling by a double cusp potential 
barrier 

V(x) = Voexp(— \x\/a) + V\exp{— \x — b\/d) (14) 

The Dirac equation (j2j) in the presence of the potential (|14|) cannot be solved exactly, 
therefore in order to discuss tunneling effects we proceed to compose the scattering 
matrices associated with the cusp at x = and x = b. This composition gives a very 
good approximation to the exact result when the separation b between the centers of 
the cups is large in comparison the shape parameters a and d. 

We proceed to calculate the scattering matrix S\ associated with the cusp potential 

A* 1 = Koexp(-M)^ (15) 

CL 

We use the following representation for the Dirac matrices [UJ 

7° = ^ 2 , l l =o-\ (16) 
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related to the matrix representation ([T]) via the unitary transformation T 

T=^=(l-za 1 ) (17) 

We have that the components of the spinor \l/ solution to the Dirac equation (T5]) in the 
presence of the cusp potential are 

— i[eVoexp(— — ) — E]) + m^f 2 = (18) 
ax a J 

^ + i[eV exp{-—)-E])^2 + m^ l = (19) 
ax a J 

where and ^2 represent the upper and lower components of ty. The solutions of the 
system of equations (}T8~|) - (fT9~}) can be expressed in terms of Whittaker functions Mk^{z) 
[18J. We solve the Dirac equation in the potential (|15p considering an incoming plane 
wave as x — >■ 00. The potential (fT5l) induces boundary conditions on the solution to the 
Dirac equation at x — 0. We split the solutions to Eqs. (fl8ft - (TT9~j) into two sectors: I 
(x < 0), and II (x > 0). Looking at the asymptotic behavior of the Whittaker function 
M k>li (z) as z -»• [IH]: 



M fciM (z) -> z ^exp{~) (20) 

we have that the solution behaving asymptotically as an incoming wave from the left 
can be written as 

* of y~ hM >*Ay) \ (9U 

Vine - a+H+k)^ (2i) 



ma 



where y = 2iaeVoexp(x / a) . The normalization constant 9 is 



m 



e = = (22) 

(2ieaV yJ2E(E + V^ 2 - m 2 ) 



The spinor \l/j nc exhibits the following asymptotic behavior as x — > —00 

1 

(i+^+fc) 



Vine ->• ( (iL +fc ) ] (^aeKofea^v^ 2 - m%) (23) 



The reflected wave is 



. 1/ *M fe _ M (y) 



where the normalization constant cj) is 

m 



(25) 



2 — m 2 



(2ieaV )-»J2E{E - VE 
the asymptotic behavior of \l/ re / as x — >• —00 is: 



#re/ ( (i-ju+fe) ] {i2aeV )~^exp{-W E 2 - m 2 x) (26) 
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The transmitted wave is : 

%rans = ( ^f ^-M ) (27) 

where y = 2iaeV exp( — x/a). The asymptotic behavior of \P trans as x — > oo takes the 
form: 

/ {\-u+k) \ 

^ trans -M "J a I (i2aeVo)"exp(i VE 2 - m?x) (28) 

which corresponds to an asymptotic plane wave traveling to the right. 

Using the spinors ^i nc ( l2~TT) . \& re / (1241) . and ^/t r (I2TI) we obtain that the solutions to 
the Dirac equation in the presence of the cusp potential in the region I (x < 0) and II 

(x > 0) are 

= M _(tet«^i^<» + M _(iz^.-iM fc+1 ,_ M ® ( 29 ) 

\ ma " /\ ma " / 

/ (|-M+fe) -|M fc+] _„fa) \ / (|+M+fc) -^M, +1; ,fa) \ 

where ai, &i, ci, di are constants. The continuity of the spinor at x = results in the 
following system 

aidM^iy) + &i<P4,-» = Cl ct> ^ +k ~^ M fc + i,- M (r/) + rfifl ^ + - + ^ M fc , M (z/) (31) 

ma ma 

(| + £; + //) (\ + k-n) 



ai9M k+ltli (u)-^ ^b l( j)M k+1 ^(v) = Cl 0M fc _ A ,(z/)+d 1 0M fc , |t (i/)(32) 

ma ma 

Expressing the outgoing amplitudes &i and d± in terms of the incoming amplitudes a\ 

and ci), we have: 

b\ = r\a\ + t[di (33) 
Ci = tiOi + r^di (34) 
where the components ri, r' 1; ti, t\ of the matrix 5i for the cusp barrier are 
r = r , = g (l+ Sa 7 ^ M ^A^)M k+ i,-M + M fc »M fc ,_» 

m J a J fe+1,— /i k,—u\ J 

m 2 a 2 JW fc+l + m k-a\ U ! 

Applying the translation property of the scattering matrix, we obtain that the scattering 
matrix S2 associated with the cusp barrier V\exp{— \x — b\/d) has the following 
components: 

^ G+"-l) 2 M 2 +M 2 (A [ ' 

m 2 d 2 J-vi K+1 _ v -+- ivi K _ r) \ K e) 
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, = _ -apbfl* {i+ S d ^ 2 M^ ;n (e)M K+l ^(e) + M K , v (e)M^ v (e) 

, _, _ 0* ^S^M^M^je) - ( -^M K , ri (e)M K+1 ^(e) 

^ (i±^l! M 2 +M 2 (e) {6J) 

m 2 d 2 m K+x _ n -t- m K _ v \t) 



where 



and 



k = iEd — ~, 77 = zciV^ 2 - m 2 , e = i2dV u (40) 



n* m i* m ,ai\ 

6 = , = =======. (41) 



{2iecV 1 ) r i^2E(E + V ' E 2 - w?) (2iedV 1 )- T iyj2E(E - VE 2 - m 2 ) 

Using the expressions ( l35l) . ( |36l) . ( |37l) . and ( |38l) we calculate the scattering matrix ( fTUl) 
for the double cusp f[T4|) . Using the components r and t of the S matrix we calculate 
the energy resonances. 

The solid line in Fig. |5] shows the behavior of the transmission coefficient T versus 
the energy for the potential f[T4"]) with a — d — 0.4 and 6 = 4. The solid line depicts de 
case when both cups have the same height Vq = V\ = 6.4271 which is a transmission 
resonance value for an energy of E res = 1.3 The solid line has four peaks corresponding 
to the scattering matrix poles E r = 1.2290 - 0.04517, E 2 = 1.3404 - 0.073917, 
E 3 = 1.79129 - 0.04237, and E A = 2.6031 - 0.03537. Fig. 1 depicts the double cusp 
configuration with Vq = V\ = 6.4271, a = d = 0.4 and 6 = 4. It is also shown the 
position of the resonances for ?RE < 3. It should be noticed that in Fig. [5] the peak 
around 1.3 does not correspond to a single resonant value but to a a superposition of 
E\, Ei and the transmission resonance E res (see Fig. H]) that modifies the standard 
Breit-Wigner profile of a energy resonance p2]. The dashed line depicts the behavior 
the transmission coefficient against against E for Vq = 6.4271 and V\ = 6. In this 
case there is not overlapping of energy resonances around E — 1.3 since there are not 
transmission resonances for cusp potentials with unequal strengths. 

Fig. 6 shows the Wigner time-delay for the two-cusp barriers described by Fig. |5l 
The solid line shows that the transmission resonance located at E res = 1.3 modifies the 
shape of the peak corresponding to the energy resonance E 2 . The dashed line shows the 
Wigner time-delay for Vq = 6.4271 and V\ = 6. 



4. Concluding remarks 

In the present article we have studied resonant tunneling of Dirac particles by a 
double square barrier and a double cusp potential. We have shown that the presence 
of transmission resonances associated with the square barrier or the cusp potential 
modifies the shape and distribution of the energy resonances. We have shown that the 
Wigner time-delay of the energy resonances is modified when they are located close to a 
transmission resonance. We have also shown that there are no transmission resonances 
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x 



Figure 4. Two cusp potentials with a = cl = 0.4, V = V\ = 6.4271 and with peaks 
separated a distance b — 4. a) The dot line corresponds to REi = 1.229, b) the 
dashed line corresponds to the transmission resonance E res = 1.3, c) The dash-dot 
line corresponds to dtE 2 = 1.340 d) The long dashed line corresponds to dtE 3 — 1.791 
e) The space dot line corresponds to 5RE4 = 1.603 Notice that E res , 3tE\ and K-Eb are 
very close each other. 

when the cusps or barriers have different strengths. The transmission with no reflection 
across a barrier whose amplitude is stronger than the energy of the tunneling particle 
is a relativistic phenomenon with no analogue in the Schrodinger framework. 

The application of semianalitical methods such as the phase- integral approximation [T9] 
or the uniform approximation [20 j could be applied in the study of transmission reso- 
nances of composite systems where no analytic solutions for the potentials are available. 
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